SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS, JUNE 2010 



1 



Adaptive Demodulation in Differentially Coherent 
Phase Systems: Design and Performance Analysis 

J. David Brown^, Jamshid Abouei^t, Member, IEEE, Konstantinos N. Plataniotis^^, Senior Member, IEEE and 

Subbarayan Pasupathy^^t^ Life Fellow, IEEE 

^ Ottawa, Canada, Tel: 613-236-1051, Email: david_jw_brown@yahoo.com 
The Edward S. Rogers Sr. Dept. of Electrical and Computer Engineering, 
University of Toronto, Toronto, Canada, Emails: {abouei, kostas, pas}@comm.utoronto.ca 



Abstract 

Adaptive Demodulation (ADM) is a newly proposed rate-adaptive system which operates without requiring 
Channel State Information (CSI) at the transmitter (unlike adaptive modulation) by using adaptive decision region 
boundaries at the receiver and encoding the data with a rateless code. This paper addresses the design and 
performance of an ADM scheme for two common differentially coherent schemes: M-DPSK (M-ary Differential 
Phase Shift Keying) and M-DAPSK (M-ary Differential Amplitude and Phase Shift Keying) operating over AWGN 
and Rayleigh fading channels. The optimal method for determining the most reliable bits for a given differential 
detection scheme is presented. In addition, simple (near-optimal) implementations are provided for recovering the 
most reliable bits from a received pair of differentially encoded symbols for systems using 16-DPSK and 16- 
DAPSK. The new receivers offer the advantages of a rate-adaptive system, without requiring CSI at the transmitter 
and a coherent phase reference at the receiver Bit error analysis for the ADM system in both cases is presented 
along with numerical results of the spectral efficiency for the rate-adaptive systems operating over a Rayleigh fading 
channel. 

- The first author completed his Ph.D. in ECE department at University of Toronto. This work was supported in part by a postgraduate 
scholarship from the Natural Sciences and Engineering Research Council of Canada (NSERC), an Industry Canada Fessenden postgraduate 
scholarship, and in part by an Ontario Research Fund (ORE) project entitled "Self-Powered Sensor Networks". 

- The material in this paper was presented in part in the 18*'* Annual IEEE Int. Symposium on Personal, Indoor and Mobile Radio 
Communications (PIMRC'07), 2007 [1]. 
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Index Terms 

Adaptive Demodulation (ADM), decision schemes. Differential Phase Shift Keying (DPSK), Differential Amplitude and 

Phase Shift Keying (DAPSK), rateless codes. 

I. Introduction 

Rate-adaptive solutions such as adaptive modulation and incremental redundancy have proven to be 
very effective at increasing the spectral efficiency of wireless systems operating over a variety of channels 
[2]-[4]. These techniques allow mobile communications systems to operate efficiently over a wide range 
of channel parameters while maintaining a target Bit Error Rate (BER), as opposed to a non-rate-adaptive 
system which is designed for a worst-case "fixed-rate" channel. In a traditional implementation of adaptive 
modulation [5], the transmitter dynamically adjusts the current level of modulation based on feedback of 
the observed Channel State Information (CSI) from the receiver. The requirement of CSI at the transmitter 
is one of the major impediments of adaptive modulation systems, especially for a high-speed mobile 
receiver for which the channel information changes rapidly. Incremental redundancy [6] systems do not 
require CSI at the transmitter, but instead encode data packets using a low-rate "mother code" and break 
the message into smaller sub-packets, sending them one-by-one until the receiver acquires enough packets 
to successfully decode the message. These systems require the receiver to decode after each sub-packet 
is received to check for successful transmission and also require a buffer large enough (in poor channel 
conditions) to accommodate the entire collection of sub-packets, meaning that the complexity of the 
decoding operation varies with the quality of the channel. 

Adaptive Demodulation (ADM), first proposed in [7] and with the expanded framework in [8], offers 
a new rate-adaptive solution that avoids the need for CSI at the transmitter (unlike adaptive modulation), 
has a fixed buffer size (roughly equal to the size of the message), and requires the receiver to perform a 
decoding operation only once (unlike incremental redundancy). In an ADM system, the transmitter sends 
data at a fixed-rate using a standard constellation and the receiver demodulates data at a non-fixed-rate 
using appropriately designed sets of decision regions chosen based on the CSI observed by the receiver. 
Before transmission, each A; -bit message is encoded using a rateless erasure code (e.g. the Luby Transform 
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(LT) code [9] or the Raptor code [10]). To recover any A;-bit message, the receiver simply buffers the 
demodulated (non-erased) bits until it has accumulated {1 + e)k "reliable" bits (where £ is a small fixed 
quantity) at which point the original message can be decoded, regardless of the erasure pattern introduced in 
the message. In [7], the ADM-based receiver uses sets of specially designed decision regions to demodulate 
some or all of the bits in a transmitted symbol (depending upon the desired level of reliability), and assumes 
a "coherent phase reference" is available at the receiver. The optimal constellations and mappings for the 
ADM system over an additive white Gaussian noise (AWGN) channel and its performance for AWGN 
were derived in [7]. Particular instances of the ADM solution are the application studies considered in 
[11] and [12] for Gaussian relay channels. 

In this paper, we consider how to implement the ADM system using differentially coherent detection: 
an attractive alternative when a reliable phase reference is not available. This is certainly a valid and 
interesting problem when the channel changes rapidly making phase coherence difficult or expensive to 
achieve. There are several works in the literature that study the performance analysis of differentially 
coherent detection in wireless RF applications (e.g., [13]) with focus on mobile satellite communications 
[14], and optical communication systems [15]. Of interest is to utilize the Differential Phase-Shift Keying 
(DPSK) in mobile communications systems which commonly circumvents the ambiguity in the phase 
recovery [16]. 

To the best of our knowledge, this is the first paper addressing this rate-adaptive technique with focus 
on differentially coherent detection in time-varying wireless systems. This paper studies the design and 
performance analysis of a hard-decision version of the ADM scheme for two common differentially 
coherent schemes: M-DPSK and M-DAPSK (M-ary Differential Amplitude and Phase Shift Keying) 
operating over AWGN and Rayleigh fading channels. The paper makes the following contributions: 

• The optimal method for determining the most reliable bits for a given differential detection scheme 
is presented. 

• Simple (near-optimal) implementations are provided for recovering the most reliable bits from a 
received pair of differentially encoded symbols for systems using 16-DPSK and 16-DAPSK. The 
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Fig. 1. Block diagram of a differential ADM transmitter/receiver. 

new receivers offer the advantages of a rate-adaptive system, without requiring CSI at the transmitter 

and a coherent phase reference at the receivers. 
• Bit error analysis for the ADM system in both cases is presented along with numerical results of the 

spectral efficiency for the rate-adaptive systems operating over a Rayleigh fading channel. 
The rest of the paper is organized as follows. Section |ll] provides an overview of the system under 
consideration and briefly introduces some ADM concepts and terminology. Section Ulllpresents the optimal 
method for determining the most reliable bits for a given differential detection scheme in the ADM-based 
system with 16-DPSK. With a similar argument as for the DPSK-ADM system model, Section HVl deals 
with a simple (near-optimal) implementation for recovering the most reliable 1, 2, 3, or 4 bits from a 
received pair of differentially encoded symbols for systems using 16-DAPSK. Sections UlI] and |IV] also 
include the Log-Likelihood Ratio (LLR) computation of each bit for every pair of received symbols. 
Section |V] includes the probability of error analysis for the near-optimal detection schemes and presents 
the spectral efficiency of ADM using differential detection. Finally, Section |VI] provides a brief summary 
and conclusions. 



IL System Model and Assumptions 

In this work, we consider an ADM-based single-hop wireless system depicted in Fig. [H where the 
transmitter desires to send data packets of equal length toward the corresponding receiver using a prede- 
termined rateless code and differential encoding. To send a fc-bit data packet, the transmitter first encodes 
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Fig. 2. Mappings of differential information for 16-DPSK and 16-DAPSK. 

it using an LT (Luby Transform) code to produce a stream of coded bits. Such codes are a class of sparse- 
graph "rateless" erasure codes which are fully explained in [9] and [10]. It is shown in [17] that LT codes 
offer strong benefits in using in practical wireless networking applications with dynamic distance and 
position nodes. LT codes have a very simple encoding mechanism that allows the transmitter to produce 
a potentially limitless sequence of encoded bits from the original k message bits. Once the receiver 
successfully recovers any (1 + £)k reliable encoded bits-where e is determined by the reliability of the 
recovered bits-a message passing decoder using the sum-product algorithm [18], can be used to recover 
the original /c-bit message. 

The coded bits are used to produce differentially encoded symbols. In this paper, we consider two 
common differentially coherent schemes: 16-DPSK and 16-DAPSK. For 16-DPSK, the phase difference 
between successive symbols conveys 4 bits of information (here denoted by 6o^i^2^3) in the standard 
fashion of [19] with the mapping shown in Fig. |2l where 6jG{0,l}, i = 0,l,2,3. Differential Amplitude 
Phase Shift Keying is a popular alternative to DPSK and has a performance advantage over DPSK at 
the expense of slightly increased constellation complexity [20]. 16-DAPSK is implemented as described 
in [20] where the amplitude difference between successive symbols conveys 1 bit of information (here 
denoted by bo) and the phase difference conveys an additional 3 bits (6162&3), according to the mapping 
in Fig. |2l By convention, a change in amplitude (i.e., consecutive symbols are transmitted on different 



SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS, JUNE 2010 



rings) represents a "1" bit, while no change in amplitude represents a "0" bit. The differential symbols are 



transmitted over an AWGN channe 



J. For the remainder of the paper, we observe the following notation. 



DPSK: After matched filtering, a sample of a 16-DPSK symbol at the receiver can be written as 

Vk = SkC^^'^ + rifc, (1) 

where = e^'^'' represents a transmitted symbol with phase 0^ and normalized symbol energy, 9^ 
represents the arbitrary (and unknown) phase introduced by the channel, and Uk is a zero-mean circularly- 
symmetric complex Gaussian random variable with variance cr^ = ^ in each dimension. It is assumed 
that 9k is constant over at least two symbol intervals such that 9^ = 9k-i. In addition, we denote the 
phase difference between consecutive 16-DPSK symbols Sk and Sk-i by 

^4>k = 4>k + 0k- (0fc-i + Ok-i) = <Pk- 4>k-i- (2) 

DAPSK: Similarly, adopting a notation akin to [20], we denote two consecutive received 16-DAPSK 
symbols (the (k — 1)"^ and k^^ symbols) as follows: 

Zk-i = dk-ie^^" + Uk-i, (3) 
Zk = akdk-ie^^" + rik, (4) 

where d^-i is a complex quantity representing the magnitude and the phase of the {k — 1)"^ symbol, 
ak is a complex quantity representing the change in magnitude and phase from the {k — 1)'*^ to the k^ 
symbol, and Uk and Uk-i are independent Gaussian noise samples as above. We note that \dk-i\ can take 



on one of only two values: Ai = ^ T+b!^ ^2 — i+r-2 » where i? = ^ is defined as the ring ratio 
in 16-DAPSK. The symbol energy is assumed to be normalized such that ^ ^ ^ = 1. The change in 
magnitude, \ak\, can be one of three values: 1/R, I, or R. Clearly, only certain combinations of \dk-i\ 
and \ak\ are allowable; for example, it would not be possible to have \dk-i\ = A2 and \ak\ = R, since 
this would require the k^^ symbol to have magnitude not equal to either Ai or A2. The following ordered 

'The Rayleigh fading case will be considered later in Section |V] For the fading channel model and as per the standard assumption in 
differential schemes, it is assumed that the complex channel gain is constant over at least two symbol intervals. 
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pairs (\dk-i\, \ak\) specify the only allowable combinations of \dk-i\ and \ak\: {Ai,l), {A2,l), {Ai,R), 
and {A2,l/R)- 

Decision Scheme: In the ADM system with differentially coherent detection, the receiver utilizes a 
/3-Decision Scheme (/3-DS) to recover the /3 most likely bit^ transmitted based on the observation of 
two consecutive differentially encoded symbols as will be described in the subsequent sections. Indeed, 
the receiver uses the instantaneous Signal-to-Noise Ratio (SNR) to determine which /3-DS to use to 
demodulate a pair of symbols in order to maintain a given BER, Pb, as a standard of reliability. For the 
above scheme, non-demodulated (unreliable) bits are erased. Demodulated (reliable) bits are buffered until 
the receiver has collected sufficient bits for decoding at which point the decoder returns the decoded fc-bit 
message. 

III. Optimal /3-Decision Schemes for 16-DPSK 

For a particular pair of differentially encoded symbols, the (3 most reliable bits in the output of the 
/3-DS should be the (3 bits with the largest magnitude Log Likelihood Ratios (LLRs). Thus, the "optimal" 
/3-DS is a device that computes the LLRs of each bit for every pair of received symbols, compares their 
magnitudes, and chooses the /3 bits with the largest LLRs. We begin with this construction and introduce 
several approximations that produce simple, near-optimal /3-DSs with easy implementations that do not 
require any LLR computations. 

Let Bij denote the set of all differential angles such that bit 6j = j, where j E {0, 1}. Thus, the 
likelihood ratio for bi for two consecutive received symbols y^. and y^-i can be written as 

Pr{yk,yk-i\bi = 0) Ea^^sb^.o Pr(l/fc, A0^) 
Myk,yk-i\bi = 1) EA</,,eB.,i Pr(z/fc,l/fc-i|A0£) 

where Pr(yfc, yfc_i| A0m) denotes the a posteriori probability of observing the pair of received symbols 
yk and yk~i given that differential angle A0m was encoded at the transmitter. Using a modified form of 
eqn. (9) in [21], it can be shown that 

Pr(y,,|/fc_i|A</)^) = Jo J , (6) 

^For instance, for the 16-DPSK or 16-DAPSK, P £ {1, 2, 3, 4}. 
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where /o(a;) is the zeroth order modified Bessel function of the first kind. It is revealed from © that only 
the Bessel term of Pr(?/fc, Z/fc-il is a function of A0m, thus allowing us to cancel the exponential 

leading terms that would appear in an expansion of (|5]). This results in Aj^ being the ratio of sums of 
Bessel terms, i.e., 

f\yk±yk-ie^^\ 

K{yk,yk-i) = — .A<A,i\ • C^) 



Note that /o(a;) is a rapidly increasing function of x, when x grows, and is approximated by /o(a;) ~ 
Thus, for moderate to large values of SNR (i.e. when is large), the sums in (|7]) can both be well 
approximated by a single dominant term, since the argument of the Bessel function will be in the steep 
part of the curve. To find these dominant terms in (|7]), we observe that 

IVk + i/fc_ie^^*'"| = VbfcP + bfc-iP + 2\yk\\yk-i\cos{(f)k - A0^), (8) 



where (pk denotes the phase difference between yk and yk-i- Clearly, ^ is maximized when the difference 
\(f)k — A0m| is minimized; thus the dominant term in the sums of Bessel functions will be the term with 
encoded differential phase Acpm ^ 5i,o closest to the phase difference (pk. Let A0j o denote the particular 
differential phase Acpm G -Bj,o such that — A0m| is minimized. Similarly, let A0j^i denote the particular 
differential phase Acpi G Bi i such that — A^^l is minimized. With these definitions and under the 

high SNR approximation, (|7]) can be written as 

h (^Vlz/feP + bfc-iP + 2|?/fe||?/fc_i|cos(0fc - A0i,o)) 
AhAyk.yk-i)^ ) ( . (9) 

-^0 i-^VlVk? + Iz/fc-iP + 2\yk\\yk~i\cos{(f)k - A0i,i)j 

To find the /3-DSs, the values of \ \ii{AbAyk,yk-i))\ must be compared for all bi over the range < 
<pk < 27r. To simplify this comparison, we define Acpi^rnin to be the value of A0j o or A0j i for which 
\(f)k - Acpi^rninl is minimized and Acpi^rnax to be the value of A0i,o or A0i,i for which |0fc - A(j)i^rnax\ is 
maximized. For example, if < 0^ < 7i"/16 then A02,mm = A02,o = and A02,max = A02,i = — 7r/4. 
Some thought reveals that all Acpi^rnin will be equal for all bi for any given value of (pk (i-C-, A0o,mm = 
A0i,mm = A02,»nm = A03,„j„). Howcvcr, tMs is not true in general for A^j^^^J^. 

^Continuing the example of < < tv/16, we will have Acfii.min = for all i, and A(j>o,max = 7r/8, A(j)i^rnax = — 7r/2, A(j)2,max = 
—n/4 (as above), and Afjjs^max ~ — tt/S. 
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Since 



^-^/bfcP + bfc-lP + 2|?/A:||yfc-l| COs(0fc - A(j)i^rnin) 
'^\J\yk? + bfc-lP + 2||/fc||l/fc-l| COS - A0i,„„^) 



> 



(10) 



then we can write 



ln(A6,(yfc,i/fc-i))| 



1 



In ( Jo ( ^y'lz/feP + Iz/fc-ip + 2\yk\\yk~i\ cos(0fc - 
1 



In /, 



'0 I 2 



yfeP + Iz/fc-ip + 2|?/fe||?/fc_i| cos I 



(11) 



Noting the fact that all A(f)i^min are equal for all i, comparing the values of | \n{Ai,^{yk,yk-i))\ can be 
achieved by comparing In (^Iq (^A/|z/fcP + |z/fc-ip + 2\yk\\yk-i\ cos (0^ - Acpi^rnax)^^ for all i. Finally, 
using the fact that ln(x) and loix) are both monotonic increasing, it can be shown that this comparison 
reduces to the following simple decision rules: 



\(j)k - A(j)i 
\(f)k - A(f)i 



> 



-Acf) 



J, max 



\\n{AbXyk,yk-i))\ > |ln {Ab^{yk,yk-^ 



< \4>k - A(j)j^rnax\ |ln (Afo^l/A;, Z/fc-i)) | < | In {Ab^{yk, yk-i))\ ■ 



(12) 



(13) 



Using the decision rules in (fT2)) and (fT3l) . it is a simple matter to determine the /9-DSs to recover the 
most likely 1, 2, and 3 bits for each pair of 16-DPSK symbols. The receiver simply computes the phase 
difference (pk between consecutive symbols and uses this angle to demodulate (3 bits based on the regions 
given in Fig. [3l where the regions are derived from the decision rules above. 



IV. Optimal /3-Decision Schemes for 16-DAPSK 

Recalling the notation from ([3]) and dH), the standard method used to demodulate 16-DAPSK is to 
compute the decision statistics = -r^^ and ipk = aig(zk) — arg(zk-i). Typically, the statistic is used 
to make a decision on the amplitude bit (bo), while the tpk statistic is used to decide on the phase bits 
(616263). We determine the optimal /3-DSs of 16-DAPSK to find the most likely 1, 2, or 3 bits based on 
the observed and ipk. 




Fig. 3. /3-decision scheme for 16-DPSK to recover 1, 2, and 3 bits for a pair of differential symbols. 

Proceeding in a similar fashion to Section HIH define Cij to be the set of all triples \ak\, arg(afc)) 

such that bit bi = j, where j E {0, 1}. The likelihood ratio for bi for the received statistics and ipk is 
thus given by 

Pr(rfc, V'fcl&i = 0) 



Fr{rk,^pk\bi = 1) 



(14) 



E(M,_,|,|a,|,argK)).eC.aP^('^^'^^l ( 1 4-1 1 , | «fc | , arg(a,) ■ 
Using eqn. (3.21) to (3.30) in [22], it can be shown that 

S P 

Pr(^fc,^fc| (|4-i|, l^fcl, arg(afc))) = / ap^^a (^^ + ^) ^fc' (^^^ 
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where 




A useful property of ([14)) when the a posteriori probabilities are given as in (fTSi) is that A;,, (rkjipk) has 
the same value for as for — . Thus, any /3-DS for 16-DAPSK can compute r', = min(rfc, — ) < 1 
and use this value for detection without loss of accuracy. This simplifies the decision schemes since the 
analysis can be confined to values of contained within the unit circle. Intuitively, it makes sense that 

Tfc and — would lead to the same reliability since it should not matter the order we transition from one 

rk 

ring to the other (since either way we will have one point on the large ring and one on the small ring). 

Note that is a rapidly increasing function of x, when x grows. Thus, for moderate to large values 
of SNR and with a similar argument as for 16-DPSK, it turns out that the sums in the numerator and 
denominator of (fT4l) can be approximated by a single dominant term. Further simplifications arise by 
noting that the denominator of (fTSi) and the multiplicative term are common to all terms in (fT4l) and 
thus will cancel. Additionally, at high SNR we note that ^ B, so we drop the additive B term 
from (fT5l) . Thus, the numerator and denominator of the likelihood ratio in (fT4l) can be approximated by 
appropriately chosen dominant terms of the form ^^es" With some thought, it is possible to identify 
the values of la^l, and arg(a;fc) that maximize (fT5l) for any particular pair of {rk.i^k)- 

A. Simple Estimate for Choosing the Differential Amplitude Threshold 

Before determining the /3-DSs for 16-DAPSK, we develop an interesting result obtained from the high 
SNR estimate of (fT4l) . In standard 16-DAPSK systems, the decision for the "differential amplitude bit" 
bo is to choose = when > Aq and to choose &o = 1 when < Aq, where Aq is an appropriately 
chosen threshold based on the ring ratio R. An optimal value of Aq can be found by numerically evaluating 
(fT4l) for i = where Aq is the value of rj^ such that A^,, (r^, ^p^) = 1, i.e., bo = and 6o = 1 are equally 
likely. The high SNR estimate of (fT4l) can be used to develop a remarkably simple (and accurate) estimate 
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for the threshold Aq. At high SNR (and for the range < r^, < 1), the dominant terms in (fT4l) are 
Pr(rfc, (^1, 1, 0)) in the numerator and Pr(rfc, ipk\ {A2, 1/-R, 0)) in the denominator, with representing 
the differential angle closest to ipk- A final approximation can be made by observing that will always 
be within ±n/8 of differential angle ^k', using this fact, we use cosi^^k — 0) ~ 1. Thus, at high SNR we 
can write 



{^)\l + rl + 2n) 



{^)\l + rl/R' + 2rk/R) 

(20) 



(a) /I + rfc^ ^ 



where (a) comes from i? = 4^, and 



To find Ao, we want to find the such that the estimate in (fT9l) is equal to 1. At high SNR, the 
behavior of the exponential determines, to a large extent, whether A{,(, (r^, i'k) is greater than or less than 
1. That is, to a reasonable approximation, ^ 1 when > Aq and ^ 1 when < Aq (at high 
SNR); due to the dominant behavior of the exponential, sign(Z) gives a reasonable estimate of whether 
or not AfcQ {r^jipk) exceeds 1. Consequently, a simple estimate for Aq is found by solving for when 
Z = 0. Solving dnj when Z = 0, results in 

A„«r,|_ = ^. (22) 

This approximation has been verified to be accurate to within two percent of the optimum value obtained 
through numerical solution [8] (over the useful range 1.5 < i? < 2.5). For R = 2, results in [23] (and 
numerical evaluation of (fT4l) ) report an optimal value of Aq = 0.68 while the estimate gives Aq ~ 2/3 = 
0.667. 



B. Optimal (3-Decision Scheme 

To compute the /3-DSs for 16-DAPSK, we first identify the dominant terms in the numerator and 
denominator of (fT4l) . Then, we use these estimates to compare the likelihood ratios of each bit in a 
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similar fashion to the way Aq was determined above. Adopting a notation similar to what was used in 
Section HIH let A^^^o denote the particular phase arg(Q;fe) (where {\dk-i\,\ak\,arg{ak)) G Q^o) such 
that — A^j ol is minimized, and define A^j i in the same manner. Furthermore, denote by Aipi^rnin 
the value of A^j o or Aipi^i for which — A^/'i,mm| is minimized, and denote by A^j maa; the value of 
or A^/'i^i for which — IS.-\\}i^max\ is maximized. 
An example computation comparing likelihood ratios is given as follows. For Aq < < 1, it can 
be shown that for % = 1,2, or 3, the numerator of (fT4l) is dominated by Pr(rfc,^fc| (Ai,l,A^jo)) 
and the denominator is dominated by Pr(rfc,^fc| (Ai, 1, A^/^j i)). Additionally, it can be shown that for 
i = 0, the numerator is dominated by Fr(rk,ipk\ (^i, ^, Aipo,min)) and the denominator is dominated by 
Pr{rk,^/Jk\ (^2, AV^cmin))- Now since one of Pr(rfc,V^fc| (Ai, 1, A^/'i,o)) and Pr(rfc, (Ai, 1, A?/'i,i)) 
must be equal to Vr{rk,ipk\ (^i, I, I^ipi^min)) (with the other equal to Pr(rfc,^fc| (Ai, 1, A^i^^^^.)))' then 
comparing \ \n{ki,^{rk.,il)k))\ with | ln(Ab^(rfc, ^/^fc))! {i ^ 0) amounts to evaluating the difference 

I ln(Afe^(rfc, Vfe))| - I ln(Abo(rfc,^/'fc))| 
= ln(Pr(rfc, V'fel (Ai, 1, AV'^,„^^„))) - In (Pr(rfc, {Ai,l, I^ipi^^ax))) 

- (In (Pr(rfc, i)k\ (Ai, 1, A^'o.m^n))) - In (Pr(rfc, ^fc] (As, l/R, AV^o,mm)))) 
= ln(Pr(rfe,^fe| (As, 1//?, AV-cm*™))) - In (Pr(rfc, (Ai, 1, A^i,^,,))) , (23) 

where (a) comes from the fact that Pr(rfc, i)k\ (Ai, 1, AV^cmm)) = Pr(r-fc, ipk\ (Ai, 1, Aipi^rnin)), and noting 
that Aipi^rnin ^rc equal for all i E {0, 1, 2, 3}. Finding the sign of the difference in (|23] ) can be computed 
by evaluating the quotient 

^ A_ Pr(rfc,y^fc| (^2, A^o,^^„)) (a) -R^ + + 2rkRcos{ipk - i^o,min) ^24) 
P^{rk,ipk\ {AiA,^i^i,max)) l + rl + 2rkCOs{iJk-'ipi,max) 

where Z = rk{Ai/ a'^Y{2R — 2cos(^/'fe — ipi^max) — T"k{R'^ — l))^ and (a) follows from the estimate 

cos{^jjk — i>i,min) ~ 1- Using the same arguments made when evaluating the quotient in (fT9l ), at high SNR 

sign(Z) gives a reasonable estimate of whether or not L > 1. This results in the following decision rule 

(when Aq < ryfc < 1): if < -^^^{R — cos{'ipk — il^i,max)) then 6j is more likely than bo; otherwise, bo is 

more likely. 
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Using this technique to compare other likelihood ratios over other ranges of we obtain the following 
simple set of decision rules to determine, for a given (rk,ipk), which bits are the most reliable: 

• Comparing \\n{AbXrk,i'k))\ with | ln(Afc^. (r^, ^fc))| /or i,j G {1,2,3}, i ^ y. 

iV'fc - AV'i,max| > l^fc - AV'j-maxI | ln(Ab^ (ffc, | > | ln(Afe . (r^, | . (25) 

• Comparing \ \\y{k],^{rk,^k))\ with | ln(Afc^ (r^, V^fc))] /or i G {1,2,3}: 



i) When > A 



0- 



rk > -I \H^bo{rk,W)\ > \HK{rk,W)\- (26) 

ii) When < Aq: 

rk > ^2 _ I ^ \ HK{rk,W)\ > I ln(Ab,(rfc, (27) 

where Aq is the decision threshold for bo. 

The decision rules in (l25l) to (l27l) are intuitively satisfying. When comparing only bits 1, 2, and 3 (i.e., 
the bits determined by the differential angle), rule (l25l) is identical to the decision rule for 16-DPSK (rule 
(fT2l)). as one might expect. When Aq < r^, < 1, rule (|26l ) indicates that as {^pk — ipi,max) increasecl (i-C-, 
bi becomes more reliable), then the right hand side of (l26l) will increase, meaning that must be large 
(i.e., close to 1) in order for bo to be more reliable than bi. Finally, when < Aq, rule (|27l) indicates 
that as {ijjk — i>i,max) increases (i.e., b^ becomes more reliable) then the right hand side of (1271 ) decreases, 
meaning that must decrease (6o becoming more reliable as gets further away from Aq) in order for 
bo to be more reliable than 6j. 

Despite the relative simplicity of the rules in (l25l) to (l27l) (compared to calculating and comparing 



LLRs), it turns out that when the rules are evaluated over all possible (rfc, il)k) they do not produce /3-DSs 
that are particularly simple to implement, unlike 16-DPSKlj. In the next subsection, we propose simple 
/3-DSs for 16-DAPSK based on a heuristic analysis; subsequent judicious application of the rules in (|25l) 
to (1271 ) allow us to lend some mathematical rigor to the simplified decision schemes. 

''An increase in {ipk — il'i,max) means that the nearest differential angle that will result in an error has become even further away from 

the observed differential angle. This indicates that the reliability of bi has increased. 

^The /3-DSs so produced have non-linear boundaries specified by the polar form equations in i26\ and jlll 
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C. Simple ^-Decision Scheme for 16-DAPSK 

Recalling from the definition of the threshold Aq, we observe that when is close in value to Aq, the 
bit 60 is very unreliable. In this case, our /3-DS (for (3 < A) should always drop ho as the most unreliable 
bit (we will specify below precisely how close needs to be to Aq for this to apply). Thus for "close" 
to Aq, decisions on 61, 62, and 63 can be made using (|25l) . resulting in a simple set of /3-DSs (similar 
to the DPSK regions in Fig. |3]). Conversely, 60 is at its most reliable when is "close" to 1 or is 
very small. So for these values of r^., the /3-DS should always demodulate 69 and demodulate (3-1 bits 
of hi, h-2, and 63, again using (|25l) . Finally, we propose a "transition region" that applies for in some 
region between 1 (i.e., 60 reliable) and Aq (i.e., unreliable) and also for in a region between Aq 
and (i.e., 60 very reliable once again). This transition region can be viewed as a hybrid of two regions: 
one region where 60 is always erased and another region where 60 is always kept. The decision scheme 
for the transition region is derived by careful consideration of the behavior of the optimum decision 
schemes that would be obtained by directly evaluating the likelihood ratios of (fT4l) over all {rk,i'k) (or 
directly implementing the rules in (l25l) to (127)) ). In actual fact, the optimum decision schemes transition 
gradually between the "60 reliable" and %q unreliable" regions in a continuous manner, passing through 
our proposed "transition region" for certain particular value(s) of r^. Intuitively the proposed transition 
region is a satisfying "halfway" point between 60 reliable and 60 unreliable. 

The simplified /3-DSs are given in Fig. |4] for /3 G {1,2,3}, where the reliable" region appears 
on the far left, the "transition region" appears in the middle, and the "60 unreliable" region appears on 
the far right for each /3-DS. Fig. |4] also introduces the threshold terms A^ 1, A^ 2, ^/3,3» and A^ 4 that 
are used by the decision schemes that recover /3 bits. For A^ 1 < < 1 or < A^ 4, bo is sufficiently 
reliable that it is always kept, while for A/3^3 < < A/3^2, &o is sufficiently unreliable that it is always 
discarded. Also, for all other values of r^, the transition region is used. A method is presented below to 
find Afs^i, A^_2, ^0,3, and A/3^4 using (|26|) and (ITTI) . 

Consider the 3-DS shown in Fig. |4l As begins decreasing (from 1), an optimal 3-DS would begin to 
transition from the %o reliable" region in a continuous fashion to the hybrid (middle) region. Fig. [5] shows 



SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS, JUNE 2010 



16 




Fig. 4. /?-Decision Scheme for 16-DAPSK to recover 1, 2, and 3 bits for a pair of differential symbols. 

this progression; Fig. [Sj^a) shows a section of the reliable region, while Fig. [5tb) shows the beginning 
of a transition to the hybrid region. In Fig. (S^b) we see that over the angles for which 63 is most reliable 
(as determined from rule (lISl)). 60 is dropped in favor of 63 when r'^ decreases by a sufficient amount. 
As continues to decrease, the angle over which 63 is selected (in favor of bo) grows as depicted in 
Fig. Sc). Eventually this angle increases to such an extent that the decision region becomes the hybrid 
(middle) region shown in Fig. Od). 

For the heuristic decision schemes depicted in Fig. |4]we opt use the "bo reliable" region when the angle 
over which 63 is selected is small (as shown in Fig. Ob)). We switch to the hybrid (middle) region when 
this angle has increased exactly halfway from zero to n/16 (where n/16 is the angle of the 63 region 
for the hybrid region as shown in Fig. Od)). This halfway point occurs as shown in Fig. [Sfc), where 
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bgOO? 



Fig. 5. Progression from bo reliable to "transition region". 
?0 




b„00? 



Fig. 6. Progression from bo unreliable to "transition region". 

the angle of the 63 region is 7r/32. The bold line in Fig. \5lc) depicts the boundary between choosing 
63 and choosing 60; thus, the bold line represents the line such that Afc„(rfc, 7r/32) = Ah.^(rk,7r /32). To 
determine the for which this equivalence occurs, we can use rule (|26l) □ by direct substitution such 



that r'. 



2{B.-COs{\l)h-ll)-i,max)) 



where il)k = 7r/32 and ip. 



3, max 



Tx/A (the correct value for ijjz^max can be 



k R2-1 

determined from FigjJtb)). This value of is precisely the threshold A3_i. 

As another example, to find the threshold A3_2, Fig. |6]proves useful as it depicts the transition from the 
"60 unreliable" region to the hybrid region. Once again, Figl6tc) shows the halfway point between these 
two regions such that the angle of the "60 unreliable" region has decreased from 47r/32 to 37r/32: exactly 
halfway to the angle 27r/32 shown in FigEd). Using rule ^ we can find = 2(fl-cosgfc-V3.n,a.)) ^ ^^^^.^ 
ipk = 37r/32 and tp3,max = 7r/4. 

In a similar fashion, the remaining thresholds for all Aj ^ can be computed. Analytical equations for 
these thresholds computed in this manner are given in Table HI The equations for the thresholds are 
modified such that A^ 1 and A^j 2 can be no larger than 1 and A^ 3 and A^ 4 must be greater than 0. This 

^Notice that we use l l26t since here we are addressing the case where > Aq as we have implicity assumed that is large. 
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modification reflects the reality that < r^, < 1. The table also evaluates the threshold equations for a 
ring ratio of i? = 2 which is the optimum ring ratio for standard 16-DAPSK in Rayleigh fading, as per 
[23]. Observe that some of the thresholds are equal to 1 or indicating that some of the regions depicted 
in Fig. m are not used. Since Fig. |4] is general and applies for any R, not all regions will be used for 
every R since the thresholds themselves are functions of R by virtue of Table H 

TABLE I 

Equations to compute thresholds Ai,j for 16-DAPSK. 



A3,2 
A3,3 
A3,4 

Aa.i 
A2.2 

A2,3 
A2,4 

Ai,i 

Al,2 

A1.3 

Al,4 



Equation 



min I 1 



2(fl-cos(7r/32-7r/4)) 



mm I i, -,-^-c°s(3^^_/32-./4)) 



ma^ (0, 2(ficos(3V3_2-V4)-l) 

max (0, 2(«£2£{ZL^az£ZiKi) ) 



min I 1 



2(ii-coB(37r/32-(-7r/4))) 



min (1, 2(«-c°s(V3_2-(-V4))) 
max (0, 2(^cos(./32-_(-V4))-i) ^ 
max (0, 2(i^cos(3./32-(-V4))-i) - 



min ( 1, 2(fl-c°s(3V8-(-V4))) 



min ( 1 



2(ii-coB(37r/16-(-7r/4))) 



max (0, 2(«cos(3./16-(-V4))-l) 



max ( 



2(flcos(37r/8-(~7r/4))-l) 



for i? = 2 



0.818 
0.745 
0.509 
0.364 
1 

0.910 
0.179 



1 

1 







As a final note, we have verified that the simplified decision schemes of Fig. |4] perform nearly as well 
as the optimum (and complex) schemes obtained by numerically evaluating (fT4l) or directly implementing 
the rules in (l25l) to (|27l) . Numerical results demonstrating this fact are provided in Section |Vl 



V. Performance of Differentially Coherent ADM 

To evaluate the performance of the ADM system for 16-DPSK and 16-DAPSK, we first examine the 
BER for each /3 decision scheme in AWGN. For 16-DAPSK, the BER can be computed by numerically 
integrating the density function in (1151) over the appropriate regions of the simplified /3-DSs as defined 
in Fig. m Additionally, we choose R = 2 which is the optimum ring ratio for 16-DAPSK in Rayleigh 
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Fig. 7. Performance of the optimum decision schemes for 16-DAPSK and "simple" decision schemes. 

fading [23]. The effect of the ring ratio in the proposed ADM system will be considered later. Fig. |7] 
illustrates the probability of error for the simplified (3 decision schemes for 16-DAPSK with R = 2. Also 
shown in Fig. |7] are BER curves (obtained through simulation) for the 16-DAPSK system using "optimal" 
(3 decision schemes that evaluate and compare the exact LLRs for every pair of received symbols. For 
[3 = 1 and /3 = 2 we observe virtually no dijference in BER performance between the optimal and 
simplified /3-DSs. For (3 = 3, the simplified scheme exhibits a loss of approximately 0.6dB at high SNR 
compared to the optimal scheme. Evaluating the BER at other values of R yield similar results (negligible 
loss for /3 = 1 and (3 = 2 and a small but measurable loss for /3 = 3 at high SNR). These results suggest 
that the simplified /3-DSs perform nearly as well as the optimum (and complex) schemes obtained by 
numerically evaluating (fT4l) . 

For 16-DPSK, the probability of symbol error for a /3-DS, Ps^p, can be computed as 




(28) 




(30) 



(29) 
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Fig. 8. BER for recovering 1, 2, 3, and 4 most likely bits for 16-DPSK and 16-DAPSK (with R = 2). 

where is the density of the phase angle between two vectors perturbed by (independent) AWGN 
samples. The integrals taken in the computation of the symbol error, above, are simply the integrals over 
the regions of received angle Acf)^ for which an error would occur (based on the /3-DSs given in Fig. [3]). 
A good approximation for integrals of the form in (l28l) to (|30l) is given in [24], eqn. (46), rewritten in a 
slightly different form below: 



Using (|28]) to (|3T]) . and making the approximation that the BER for a /3-DS is ^Ps,(3, the probability 
of error curves for 16-DPSK are given in Fig. [8l The figure also provides the BER curves (from Fig. |7]) 
for 16-DAPSK with R = 2. Note the initially surprising result that although 16-DAPSK outperforms 
16-DPSK when recovering all four bits (as is well known), when recovering 1, 2, or 3 bits 16-DPSK 
shows a clear advantage over 16-DAPSK. 

The fact that 16-DPSK outperforms 16-DAPSK for lower DS can be explained by noting that in a 16- 
DPSK scheme, certain bits are very well protected while others are more susceptible to error. By erasing 




(31) 



where erfc(a;) = e dt is the complementary error function. 
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the error prone bits (for a given received angle) large gains are observed; however, for 16-DAPSK the bits 
are more uniformly reliable/unreliable with the consequence that erasing the least reliable among them 
will not result in gains as significant as for 16-DPSK. More specifically, for 16-DPSK, the least reliable 
bit (in general) is the bit that changes value between two differential angles denoting two "adjacent" 
symbols. When this bit is dropped (e.g., for 3-DS), the performance of the remaining bits is reminiscent 
of an 8-DPSK system — i.e., a significant gain is observed. With 16-DAPSK, an advantage is gained for 
4-DS since the minimum differential angle separating differential symbols is 7r/4 (as opposed to tt/S for 
16-DPSK). Of course, this increase in minimum angle is achieved by allowing for two distinct amplitudes: 
Ai and A2. This strategy does indeed succeed in reducing the error rate of the least reliable bit. Often, 
however, the least reliable bit for a 16-DAPSK constellation is in fact the differential amplitude bit; thus, 
when this is neglected, the remaining bits are less reliable than for a 16-DPSK 3-DS since the inner ring 
necessarily has a smaller radius than the standard 16-DPSK radius (assuming equal symbol energies for 
the two constellations). 

With the BER computed for all /3-DSs, it is now possible to compute the optimal operating regions for 
the ADM system over a Rayleigh fading channel (i.e., compute which (3-DS to use for a given observed 
instantaneous SNR). These operating regions are computed using the techniques outlined in [25] assuming 
that the fading is constant over at least two symbol intervals. Using the computed operating regions, the 
spectral efficiency (for "raw" uncoded bits) for ADM using 16-DPSK and 16-DAPSK is given in Fig. |9l 
It is seen that ADM using 16-DPSK outperforms ADM using 16-DAPSK for low rate transmission (until 
about 2.5 bits per symbol); additionally, it does not suffer a large performance deficit at higher rates. This 
is explained by the relatively poorer performance the 1-DS, 2-DS, and 3-DS of 16-DAPSK compared to 
16-DPSK as noted in Fig. E 

Although [23] demonstrated that a ring ratio of i? ^ 2 is an optimal choice for 16-DAPSK operating 
in a Rayleigh fading environment, it is not obvious that this ring ratio remains optimal for an ADM 
system using 16-DAPSK. Reference [8] demonstrates the effect of the ring ratio on the BER performance 
(in AWGN) for each (3 decision scheme. It is shown in [8] that for the 4-DS (i.e., standard 16-DAPSK 
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Fig. 9. Spectral Efficiency for 16-DPSK and 16-DAPSK (with 7? = 2). 

differentially coherent demodulation), the BER is minimized for a ring ratio of = 2, a result that agrees 
with [23]. However, for the 3-DS, 2-DS, and 1-DS, it observed that the BER tends to increase with R, 
exhibiting no local minimum at i? = 2 (e.g., see Fig. [10] for 2-DS). Thus, for all /3-DS except 4-DS, 
the BER performance improves for lower ring ratios. These results can be explained by observing that 
as R decreases below 2, the inner and outer ring radii become closer in value; consequently 6o becomes 
increasingly unreliable. This results in an increase in BER for 4-DS. For all other /3-DS, however, the 
unreliable 6o is most likely the bit that will be discarded; the remaining bits become more reliable for 
smaller R since the radius of the inner ring increases, leading to more reliable angle measurements. 

Since the BER of all /3-DS (including 4-DS) increases when i? > 2, we conclude that it is never 
advantageous (at least in terms of BER performance) to use ring ratios larger than 2 in the ADM system. 
The question still remains, however, what is the best ring ratio {R < 2) for ADM. The answer depends 
on the desired system performance and range of operation. For example, if the system is expected to 
operate primarily in high average SNR at close to 4 bits per symbol, then i? = 2 is a good choice since 
the system will spend most of its time using 4 — DS. However, if the system is expected to have a very 
broad range of operation, a smaller R may be advisable. 
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Fig. 10. BER of 16-DAPSK 2-DS for various ring ratios. 

In Fig. [m the spectral efficiency curves (for Rayleigh fading) are computed for a range of ring ratios. 
As expected, the ADM systems using low values of R outperform the ADM system using i? = 2 at lower 
rates; however, these systems are less effective at higher rates when the performance of the 4-DS becomes 
more significant. 

VI. Conclusion 

In this paper, we have derived optimum and near-optimum receivers for differentially coherent reception 
in Adaptive Demodulation systems using 16-DPSK and 16-DAPSK. Simple decision rules were derived 
for the optimal demodulation of the (3 most likely bits for both 16-DPSK and 16-DAPSK. Although these 
simple rules led to simple decision schemes for 16-DPSK, the resulting decision schemes for 16-DAPSK 
contained nonlinear boundaries that were not conducive to a simple implementation. Based on the 16- 
DAPSK decision rules, heuristic near-optimal decision schemes were proposed for 16-DAPSK with a 
simple implementation; in addition, formulas were derived to tailor the decision schemes for any desired 
ring ratio. The probability of error for 16-DPSK and 16-DAPSK receivers demodulating the most reliable 
1, 2, 3, and 4 bits per symbol was computed along with the spectral efficiency of the ADM systems. 
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Fig. 11. Spectral Efficiency of 16-DAPSK for various ring ratios (raw BER of 10 ''). 

A surprising result was demonstrated that over a large operating region, 16-DPSK actually outperforms 
16-DAPSK for Adaptive Demodulation systems. The impact of the ring ratio on the spectral efficiency 
of ADM using 16-DAPSK was investigated and it was shown that ring ratios larger than 2 are never 
beneficial (in terms of spectral efficiency) for such systems. A tradeoff was shown to exist where low 
ring ratios improved the performance of 16-DAPSK for low average SNR operating regions, while larger 
ring ratios led to an improvement in the performance of high average SNR operating regions. 
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